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We investigate the interiors of 3+1 dimensional asymptotically flat charged and rotating black 
holes as described by observers who fall into the black holes at late times, long after any perturbations 
of the exterior region have decayed. In the strict limit of late infall times, the initial experiences of 
such observers are precisely described by the region of the limiting stationary geometry to the past 
of its inner horizon. However, we argue that late infall-time observers encounter a null Shockwave 
at the location of the would-be outgoing inner horizon. In particular, for spherically symmetric 
black hole spacetimes we demonstrate that freely-falling observers experience a metric discontinuity 
across this shock, that is, a gravitational shock-wave. Furthermore, the magnitude of this shock is 
at least of order unity. A similar phenomenon of metric discontinuity appears to take place at the 
inner horizon of a generically-perturbed spinning black hole. We compare the properties of this null 
shockwave singularity with those of the null weak singularity that forms at the Cauchy horizon. 
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I. INTRODUCTION 



In Einstcin-Hilbert gravity coupled to various matter fields, the exterior geometry of a 3+1 dimensional asymp- 
totically flat black hole (BH) spacetime typically approaches a stationary solution at late times. Non-stationary 
perturbations decay both by falling across the horizon and dispersing to infinity, as described by the ringdown of 
quasi-normal modes followed by power-law tails. Our purpose here is to explore a corresponding late-time limit of the 
associated black hole interiors. We will argue that as far as the observations of late-infalling physical observers are 
considered, the result is well-described by a simple effective geometry which contains the part of the corresponding 
stationary BH solution to the past of the inner horizon. However, the regular inner horizon is replaced by singular 
components of two different types: (i) The ingoing section of the inner horizon — the Cauchy horizon (CH) — is replaced 
by a null, weak, curvature singularity, and (ii) the outgoing section of the inner horizon is replaced by an outgoing 
shock-wave singularity. The presence of a null, weak, curvature singularity at the CH is a well-known phenomenon 
since the pioneering works of Hiscock on the Reissner-Nordstrom-Vaidya solution [1] and of Poisson and Israel on the 
mass-inflation model [5] (see also [3]). It is the second singular component — the outgoing shock- wave singularity — 
which will be our main concern in this paper. Our study is motivated in part by the picture of extreme black holes at 
late times suggested in [5j and explored further in [5] . This picture agrees with the extreme limit of our results below. 

The starting point for our analysis is the large body of literature studying perturbations of the Reissner-Nordstrom 
(RN) and Kerr interiors. With the assumption of spherical symmetry, these works establish that perturbations 
transform at least the initial part of the ingoing inner horizon (the CH) of RN into a null curvature singularity often 
called the mass-inflation singularity. This singularity is weak in the sense that the metric remains continuous at 
the singularity, though it is not diffcrentiable. In particular, the area-radius r of the spheres is well-defined at this 
singularity, taking the value r_ = M — \J M 2 — q 2 near the point marked i + on the conformal diagram shown in figure 
[I] (left) and shrinking toward the future as described by the Raychaudhuri equation, eventually reaching r — 0. In the 
spherical case, at least when the matter content includes a minimally-coupled massless scalar field, it was numerically 
established [SHI] that when r shrinks to zero the weak null singularity meets a strong spacelike singularity 1 , along 
which r — 0. This situation is depicted in the left panel of figure [TJ Note that curvature scalars diverge at both the 
above spacelike and null singularities. While the establishment of a spacelike singularity is less certain for other forms 
of matter (e.g. perfect fluids), it is nevertheless expected that all future-directed timelike and null curves inside the 
black hole will be incomplete due to reaching a curvature singularity of some form. 

Investigations of nonlinearly-perturbed spinning BHs reveal a similar scenario. Perturbative analyses [9141 1 j again 
indicate the formation of a null weak scalar-curvature singularity at the CH (though this time the singularity is 
generically oscillatory [IT], as opposed to the monotonic mass- inflation singularity in the spherical case). The presence 
of such a null singularity is also supported by an asymptotic local analysis of the Einstein equations [12] , as well as by 
exact analytical constructions of locally-generic classes of null weak singularity [T3]. Yet, to the best of our knowledge, 
no numerical verification of this scenario has yet been carried out in the spinning case. 

The results above provide a good (if not yet fully complete) understanding of generic (charged, rotating, non- 
stationary) black hole interiors. However, the detailed experiences of any given observer inside such a black hole will 
in general depend on the process by which the black hole was formed and on the particular perturbations generated. 
In contrast, we argue below that the spacetime effectively simplifies from the point of view of observers who enter the 
black hole at late times, which we call late infall-time (or just late-infall) observers. The simplified spacetime may be 
described by the simple, stationary BH solution up to the inner horizon — and an effective outgoing shock wave at the 
outgoing portion of the latter (plus a null weak singularity at the CH). This structure is depicted in the right panel 



Due to the key role it plays in our analysis, it is useful to describe the notion of late-infall observers in more 
detail. We start with a simple demonstration (though not necessarily the most precise or most general one) of this 



These numerical simulations showed that a transition from a null to a spacelike singularity occurs in a region which (when mapped 
to a collapsing-shell spacetime) would correspond to being outside the shell (i.e. where the electric field is non-zero). However, when 
the initial scalar perturbations are sufficiently weak, the full focussing of the CH to r = will only occur inside the shell. We strongly 
expect the formation of a spacelike singularity in this case as well, and have drawn this scenario in figure [T] 



of figure [JJ 
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FIG. 1: Left: A Reissner-Nordstrom black hole formed by collapse of a thin spherical charged shell and subject to spherical 
perturbations by a massless scalar field. Two types of singularities form: A null weak singularity at the CH (the diagonal 
dashed line) and a spacelike r — singularity (the solid horizontal line; though see footnote [TJ. Depending on the strength of 
the initial perturbation, the spacelike piece of the singularity might also form earlier so that it intersects the worldline describing 
the shell. Right: Our proposed effective geometry of the perturbed (but spherically symmetric) collapsing shell spacetime as 
seen by late-time observers. It consists precisely of the region of the unperturbed eternal Reissner-Nordstrom black hole with 
r > r_, together with two types of singularities on its future boundary (at r = r_): The diagonal dashed line at the upper 
right represents the mass-inflation singularity, a null weak singularity located at the CH. The diagonal solid line at the upper 
left represents a null shock-wave singularity, where the metric tensor is effectively discontinuous. 



concept. Recall that perturbations outside the BH decay at late time. Therefore, for observers who fall into the BH 
at sufficiently late time, the exterior geometry will be well approximated by the stationary (and axially-symmetric) 
BH metric. This allows one to associate specific values of energy E (= —u t ) and angular momentum L (= u v ) to 
the late-time geodesies. More importantly, owing to the approximate time-translation of the external geometry, from 
any "seed" infalling geodesic To, we may construct a one-parameter family Y of similar geodesies, obtained from To 
by time translation to the future. (We emphasize that in the present construction the members of Y are exactly 
geodesies, all related to T$ by the approximate time translation. 2 ) Note in particular that all geodesies in Y share 
(approximately) the same values of E and L. Now, each member of Y is characterized by the parameter v e h, namely 
the value of the advanced time v (Eddington's advanced null coordinate) at which the geodesic crosses the event 
horizon. For any given seed geodesic r , the late-infall observers are those members of Y characterized by sufficiently 
large values of v e h- 

This concept of late-infall observers may be generalized, and reformulated in a somewhat more precise manner 
(though we will not attempt a fully precise definition here) . Consider a continuous one-parameter family of inextendible 
causal curves labeled by the advanced time v e h at which they cross the event horizon, with v e h taking values in some 
range (vo, +oo) so that the family includes curves that enter the black hole at arbitrarily late times. Note that since 
perturbations outside the black hole decay, the advanced time v can indeed be used as a coordinate along the horizon 
at sufficiently late times. We require that, in the limit v e h — > oo, the part of these curves to the past of the event 
horizon approaches a stationary family of such curves in some stationary spacetime; i.e. for which curves with different 
values of v e h are related by the corresponding time-translation. We further require that, for large v e h, the parts of 
our curves to the future of the event horizon all have the same proper accelerations when expressed using a reference 
frame parallel-propagated along the world line in terms of the proper time along the worldline after crossing the event 
horizon. One may think of this as the assumption that all observers in a given family are equipped with identically 
pre-programmed rocket ships. Of course we insist that these reference frames at the event horizon are also related by 
an approximate time translation at large v e h- For our rather qualitative purposes below it will not be necessary to 



2 The geodesies in T may be constructed as follows. Let x M = (t, x l ) be a set of coordinates for the BH exterior, such that at late time the 
associated metric functions are approximately independent of t. We pick a certain point Pq on the "seed" geodesic To outside the BH, 
and time-translate it to the future by a certain amount At. At the new point, which we denote Pq, we set the four- velocity components 
u M to be numerically the same as those of To at Pq. The geodesic F' which emanates from Pq with those initial four- velocity components 
now becomes a member of the set I\ 
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specify the precise rate at which these limiting behaviors are approached, though some such specification will certainly 
be needed to derive more precise results. 

After a brief review of charged spherical black holes in section [IT| w e study the experiences of late-time observers 
in stationary spacetimes subject to linear perturbations in sectiormll] Such observers experience no perturbation at 
all until they would expect to encounter an inner horizon. However, they effectively encounter a Shockwave at the 
outgoing inner horizon. While we include a discussion of linearly perturbed Kerr black holes, our treatment mainly 
focuses on the simpler spherically symmetric case. 

Non-linear perturbations are addressed in sections |IV| and |V| Here we consider only spherical black holes. Section 



IV addresses the model of a charged BH perturbed by a self-gravitating scalar field. We show that the experiences 



of freely-falling late-time observers again agree with those in unperturbed Reissner-Nordstrom up to the point where 
they would expect to encounter the (outgoing section of the) inner horizon. However, instead of finding a smooth 
null surface at that point, they effectively encounter a gravitational Shockwave at which the metric is discontinuous. 
Section [V] then gives a heuristic argument that the experiences of more general late-time observers are similar, and in 
particular that they are described by the effective spacetime shown in the right panel of figure [T] The final discussion 
in section |Vlj describes possible generalizations to rotating black holes and to black holes in any dimensions. We also 
discuss several aspects of the gravitational shock-wave phenomenon which takes place at the inner horizon. 



II. PRELIMINARIES: SPHERICAL CHARGED BLACK HOLE 



The RN solution is the unique spherically-symmetric electrovac geometry. In Schwarzschild coordinates (t, r, 9, tp) 
it takes the form 

ds 2 = -Fdt 2 + F~ 1 dr 2 +r 2 dQl, (2.1) 

where F = 1 — 2M/r + q 2 /r 2 , and dil 2 = dd 2 + sin 2 (6)dcp 2 is the unit two-sphere. Throughout this paper we shall 
consider the non-extreme black-hole (BH) case, namely < \q\ < M. In this case F(r) vanishes at two r values, 
r± = M ± \J M 2 — q 2 . The larger root r + corresponds to the event horizon, and the smaller one r_ to the inner 
horizon. Figure [2] (left) depicts a part of the Penrose diagram of the eternal, analytically-extended, RN geometry. 
Note that the inner horizon has two separate portions - the two intersecting null lines denoted "r_". 




FIG. 2: Left: The eternal Reissner-Nordstrom black hole. The vertical lines are singularities at r = 0, while the diagonal lines 
describe copies of asymptotically flat future or past null infinities (I + and outer horizons (r+), and inner horizons (r_). 
The full analytically extended solution consists of a periodic vertical array of copies of the regions shown. Right: A conformal 
diagram for the spacetime of a collapsing (but otherwise unperturbed) spherically symmetric massive charged thin shell. The 
r = curve is a regular origin in the flat region inside the shell. The Cauchy horizon is also indicated. 



Later we shall also consider the spherically-symmetric spacetime of a charged collapsing thin shell. The geometry 
outside the shell is described by the Reissner-Nordstrom (RN) metric (2.1). Inside the shell the geometry is flat, 
i.e. Minkowski. Figure [2] (right) displays this hybrid spacetime. In this non-eternal BH spacetime the ingoing 
portion of the inner horizon is a Cauchy horizon (CH) . The diagram only displays the globally- hyperbolic piece of the 
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spacetime — namely, the region up to the Cauchy horizon (the extension beyond the CH will not concern us in this 
paper). 

Let us examine the free-fall orbits of observers who jump into the BH. Throughout this paper we assume, as 
usual (and without loss of generality) that the motion is confined to the equatorial plane. These geodesic orbits are 
characterized by two constants of motion: the "Energy" E = —u t > 0, and the angular momentum L = u v . They 
satisfy the radial equation 3 

J = -y/E*-(l + L 2 /r 2 )F(r). (2.2) 
dr 

Of particular importance will be the behavior of these worldlines in the neighborhood of r — r_ , where the last result 
becomes 

In the analysis below it will be useful to express the RN metric in double-null coordinates. A particularly useful 
form is given by the Eddington coordinates 

u = t — , v = t + r* , (2.4) 

where is the tortoise coordinate defined by 

dr/dr* =F(r). (2.5) 

The line element then becomes 

ds 2 = -F(r)dudv + r(u, vfd^l , (2.6) 



where r is to be regarded as a function of u and v, determined (implicitly) by setting = (v — u)/2 in Eq. ( |2.5[ ). 
Notice that in the region which will mostly concern us here — the domain r_ < r < r + — u is past-directed (this choice 
of sign simplifies many of the expressions below). 

Note that — ¥ +oo at r — r_, implying that either u or v must diverge there. It follows that at the ingoing section 
of the inner horizon v — ¥ oo (as u is regular), whereas at the outgoing section u —¥ —oo (and v is regular). 

Since F vanishes at the inner horizon, in its neighborhood we m ay a pproximate F(r) = — 2n(r — r_), where 



K = — (l/2)(dF/dr) r=r _ (note that n > 0). It then follows 4 from Eq. (2.5 1 that near the inner horizon 

r - r_ = Me- 2Kr * . (2.7) 
That is, diverges logarithmically (in r — r_) at the inner horizon. 



The metric (2.6) is singular at the inner horizon (where det(g) vanishes). To remove this coordinate singularity 
we define the inner-horizon's Kruskal-likc coordinates U = —e KU , V = —e~ KV . With this choice of signs, both U and 
V are future-directed. Note that V (U) vanishes at the ingoing (outgoing) section of the inner horizon. The line 
element now reads 2gjjyd\J dV + r 2 d£l?,. We shall not need here the specific form of the metric functions r(U, V) and 
9uv{U, V) = guv( r )- We shall just note that both functions are regular. Furthermore, guv turns out to be a smooth 
function of r which (unlike F) is nonvanishing at r_. As a consequence, the Kruskal metric is regular at the inner 
horizon. 



III. SIMPLE EXAMPLES OF LATE-TIME PERTURBATIONS: LINEARIZED FIELDS 

In this section we consider several types of linear perturbations, and explore how these perturbations are experienced 
by late-time infalling observers. 

We shall start our analysis by addressing the simplest type of perturbation, namely a purely-outgoing, spherically- 
symmetric, test scalar field. Then we shall proceed to consider more realistic types of linear perturbations, deferring 
non-linear perturbations to the following sections. 



At a certain value of r, r = rtounce < f—, dr/dr flips its sign (a phenomenon known as "gravitational bounce"). In this paper, however, 
we f ocus on the behavior of orbits up to their first approach to r = r_, hence dr/dr < 0. 

Eq. ]2.5\ d efines r„ up to an integration constant. We use this freedom and choose the convenient pre-factor M in the right-hand side 
of Eq. (12. 7ft , which fixes this arbitrary constant. 
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A. Simplest example: monotonic, outgoing, test scalar perturbations 

Consider a free, massless, minimally-coupled, test scalar field 4> on the RN background. The general behavior of 
a scalar field of this type inside the BH will be discussed in the next subsection. Here we focus on a spherically- 
symmetric scalar held, restricting our attention to the very neighborhood of the Cauchy horizon. As it turns out, in 
this region the held becomes purely outgoing, namely 5 <f> = P(u). In our hrst example we shall further assume, for 
simplicity, that P(u) vanishes until a certain value u\\ then it grows monotonically up to u = U2, where it reaches 
its hnal value Pt (and remains constant afterward). While this type of function P(u) is certainly oversimplified, it 
transparently demonstrates the mechanism responsible for the shock-wave formation. 6 Both null lines u = u\ t 2 are 
assumed to reside in the internal range r_ < r < r+. We shall now explore the behavior of <fr along the worldline of 
the infalling observer, as a function of proper time r — focussing our attention on late-time observers. 

A free-falling observer in a RN spacetime crosses the inner horizon r = r_ through its outgoing section (where v is 
finite) . Let us denote the values of v as the observer crosses the event and inner horizons by v e h and Vih , respectively. 
Since we are interested in late-time observers we shall assume v e h S> M. 7 Along a timelike worldline v increases 
monotonically, hence in the range between the two horizons v e h < v < v^- 

Obviously, u e /, and vn L depend on the infalling geodesic. However, for hxed parameters E and L, the difference 
Av = — Veh — Av(E, L) will be independent of the infall time, owing to the time-translation symmetry of the RN 
metric. Av is typically of order M (we assume here that E is of order unity, i.e. not too large or too close to zero, 
and \L\ is ~ M or smaller; and similarly the BH is not too close to Schwarzschild or to extremality.) 

Next we evaluate the proper times T\ 2 at the two events where the worldline intersects the null lines u = 1*1,2- 
For convenience we set r = at the worldline's intersection with the inner horizon. From the above assumptions (in 
particular v e h S> M) it immediately follows that i\ M — and hence r = r_ — throughout the portion u\ > u > u 2 



of the worldline. We can therefore use Eqs. (2.3) and (2.7) in this domain to obtain 



r = -E~ l {r - r_) ~ -JfrV*"'. (3.1) 

Substituting = (v — u)/2 we hnd for t\ and r 2 

7-1,2 = -ME- 1 e K{ui - 2 - Vl - 2 \ (3.2) 

where v\^ respectively denote the value of v at the worldline's intersections with u = u\ t 2- 
Since V\ i > v e h, we readily hnd 

|ri, 2 | < (ME- 1 e KUl - 2 )e- KV " h > (3.3) 

hence the difference Ar 12 = r 2 — t\ > is bounded by 

A712 < (ir^""^-^, (3.4) 

The last inequality tells us at once that the late-time infalling observers will see a scalar-held profile 4>{t) which 
rises from zero (at «i) to its maximal value Pf > (at u 2 ) within an arbitrarily short proper-time interval, oc 
exp(— Kv e h)- For a sufficiently large v e h, this proper-time interval will presumably be unresolved by an observer with 
hxed sophistication. The large-We^ observers will thus experience the scalar perturbation as a sharp shock-wave of 
hnite amplitude. Note that in terms of the proper time of these observers, the shock wave is detected effectively at 



t = [the limit v e h — > 00 of Eq. (3.3)] — namely, just at the crossing time of the (outgoing portion of the) inner 
horizon. 



late infalling observers. Based on Eq. (3.2) we express Ari 2 as 



It may be instructive to complement Eq. (3.4) by an actual estimate of (rather than just a bound on) Ati 2 , for 



Ari2 — — (e K { Ul+ ( Vih ~ Vl ^ — e K ^ U2 +( Vih ~ V2 ')^ e ~ KVih . (3.5) 



5 This simple form follows from the behavior of perturbation fields inside BHs, which we discus s in t he next subsection: In the very 
neighborhood of the inner horizon, the outgoing and ingoing modes effectively decouple [see Eq.| |3,ld| l]. Furthermore, the u-dependent 
component of <f> decays as v~ n at large v, leaving us with only the ^-dependent component. 

6 Essentially the same analysis will apply in the more general (and more realistic) case, in which the outgoing field starts right after the 
event horizon, and may also continue its growth beyond the (outgoing portion of the) inner horizon. Also the function P(u) needs not 
be monotonic (see discussion below). Here we picked a monotonic function P(u) only for the sake of conceptual simplicity. 

7 Within the shell-collapse model we further assume that v e h is sufficiently large that the orbit is confined to the shell's exterior throughout 
the range r > r_ . Also, in the case iti > we further demand u e jj 2> M + U\, such that r, > M as m = U\ is approached. 
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Consider now the quantity Vih — v\. For fixed parameters E and L, dv/dr is a well-defined function of r (independent 
of v e h)- In the late-infall limit, r — > r_ when the orbit reaches u\. The quantity dv/dr at u\ similarly approaches 
its inner-horizon value, which (for E > 0) turns out to be a definite finite number, (1 + L 2 /r 2 _)/(2E). In addition, 
from Eq. (3.3 1 it follows that t\ — > in the late-infall limit, and thus that Vih — V\ — > 0. Obviously exactly the same 
argument applies to v^ — V2 as well. Using — «i,2 — > and = v e h 



Av in Eq. (3.5), we thus obtain 



Ar 12 = 



M , 

~e' 



")e 



-reAt 



(3.6) 



(Note that the pre- factor in squared brackets is independent of the infall time v e h-) 

The typical behavior of <j> as a function of the observer's proper-time is depicted in Fig. [3] for several values of infall 
time v eh . 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 



FIG. 3: The typical behavior of cj> as a function of the observer's proper time is depicted for several values of v e h- 



Summarizing, the late-time infalling observers will experience the scalar perturbation as a sharp shock wave of 
finite amplitude (A<fi = Pf), with effectively vanishing rising time (oc e~ KVeh ) — located just at the outgoing section of 
the inner horizon. 

So far we have considered the case of a (near-CH) outgoing field cf> = P(u) which varies monotonically from P = 
at u — u\ until u — u^ where it saturates at P — Pf . The extension of the argument to a more general function P(u) 
(not necessarily monotonic; and not necessarily one with well-defined final saturation value) is straightforward: It is 
sufficient to assume that in the near-CH region, P(u) varies over a certain range u\ > u > U2, from Pi at u\ to some 
P2 7^ Pi at U2- Then, the above analysis yields that for a late-infall observer the finite jump A<f> = P2 — Pi 7^ in the 
value of <p w iU take place within an extremely small proper-time interval, which decreases as e~ KV " h . For an observer 
of given sophistication, this exponentially small proper time will become unresolvable at sufficiently large v e h- Hence, 
a late-infall observer will again experience a kind of effective shock-wave phenomenon: a finite change in the field, 
within an effectively-vanishing proper-time interval. 8 



B. More realistic linear perturbations 



We turn now to consider more realistic perturbations of a spherical charged BH. We still assume that in the region 
of interest the perturbations are small and can therefore be treated linearly. Again, we would like to explore how 
these perturbations will be experienced by late-time infalling observers. 



1. Non-spherical test scalar field 



The linearized perturbations — both outside and inside the BH — may conveniently be toy-modeled by a minimally- 
coupled massless scalar field <fi, satisfying g a ^4>;ap — [H]- The perturbation field may be decomposed in spherical 
harmonics Y\ m in the usual way, 



Y, Y im^^)4>im{r,t). (3.7) 



A "classic" shock wave typically contains the following three components: (i) a steady pre-shock value, (ii) a (different) steady post-shock 
value, and (iii) a sharp transition between the two, through a transition region of arbitrarily-small width. Here we use a somewhat 
generalized notion of shock wave, because we do have the component (iii), and effectively also (i), but we don't necessarily have (ii). 



The individual perturbation modes cj>i m all satisfy a hyperbolic partial differential equation of the form 

-^5— -^- = V|(rM m> (3.8) 
with i\)i m = r<j)i m , and with an /-dependent effective potential given by 



V t (r)=F 



1(1 + 1) 2{M-q 2 /r) 



(3.9) 



One then finds that outside the BH all modes <j>i m decay to zero at late time [14. This decay typically proceeds 
in two stages: First is the stage of "quasi- normal ringing" (i.e. exponentially-damped oscillations). Subsequently, 
after the ringing has damped, the late-time perturbations are dominated by inverse-power tails. All modes <f>i m decay 
(along worldlines of constant r) as t~ n , with n = 21 + 3 (or n — 21 + 2 for initially-static multipoles) [13 \W\ — the same 
inverse-power form as in the Schawarzschild case [14]. Throughout the rest of the paper we shall focus our attention 
on the inverse-power tails (rather than the quasi-normal ringing), because it is this component which eventually 
dominates at late time. 

Investigations of the dynamics of linearized perturbations inside the BH reveal a behavior which parallels the 
external dynamics in many respects, though there also are some important differences. The infalling power-law tails 
lead to a similar inverse-power decay inside the BH. Thus, along lines of constant r (between the event and inner 
horizons), the perturbations still decay as t~ n [T7j. Note, however, that this time t is a spacelike coordinate, so what 
we face here is a spatial rather than temporal decay. (For a discussion of this issue see Ref. [18].) 

Being proportional to F, the effective potential Vi(r) vanishes at r_ (exponentially in r*, as it also does at r + ). As 
a consequence, in the neighborhood of the inner horizon the perturbations take the simple form of a superposition of 
outgoing and ingoing modes, that is, 

4>i m = P u {u)+P v (v). (3.10) 

Hereafter we omit the sub-indices l,m from the P-coefficients for brevity. 

In the asymptotic region h,d> M, both functions P u and P v admit simple inverse-power forms: 

P u (u)=au- n , P v (v)=bv~ n , (3.11) 

again with n = 21 + 2 or n = 21 + 3 (depending on the presence or absence of initial static multipole) . The coefficients 
a, b are determined by the scattering problem, and are given in Ref. [T7] (see also [T5]1. 

We consider now a late-infall observer, and examine how this observer will record the scalar perturbations, focusing 
on a particular multipole 4>i m - We restrict our attention to the orbit's section between two hypersurfaces u = U\ and 
u = u 2 (both reside in the domain between the event and inner horizons, and satisfy u 2 < u\). We focus on the 
quantity A<pi m = 0; m (r 2 ) — 4>i m { T i)i where as before, respectively denote the proper times when the observer 
crosses the two hypersurfaces u = iti,2- 

For a sufficiently late infall time, the observer will reach u = u\ only when the orbit is already in the near-CH 



region, where Eq. (3.101 applies. Therefore, 

A0, m £* [P u {u 2 ) - P u ( Ul )} + [P v (v 2 ) - P>i)], 

where v± } 2 = v{jx,2)- Furthermore, since P v cx v~ n at large v (and v 2 > v± > v e h), it follows that the second term in 
squared brackets vanishes at the late-infall limit. Thus, for a sufficiently large v e h the expression for A<j>i m simplifies 
to 

A^ ro £* P a {u 2 ) - P u (ui). (3.12) 

For generic choice of u\, u 2 the RHS is non- vanishing. (This is most explicitly verified in the case where U1.2 are both 
in th e ea rly domain, u 3> M, where P u cx u~ n . 9 ) Yet, the proper-time difference AT12 = t 2 — t\ is still given by 



Eq. (3.6), namely, AT12 cx e KV " h . Thus, for sufficiently late infall, the observer will watch a finite jump A(f>i m in the 



field component <^ m , during an effectively vanishing (i.e. physically non-resolvable) proper-time interval. 



Yet, a more appreciable A0; m should be achieved when U2 is not S> M (and u± is not too close to M2). 
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2. Linear electromagnetic and gravitational perturbations 

In the RN spacetime, owing to the presence of electric field, the (electrovac) gravitational and electromagnetic 
perturbations are mutually coupled already at the linear level. Still, one can write decoupled field equations for a 
pair of combined electromagnetic/gravitational variables (i.e. two specific linear combinations of the electromagnetic 
and gravitational variables) 20J. In particular, based on a formalism developed by Moncrief [3T], Gursel et al. [S3] 
constructed such a pair of electromagnetic/gravitational field variables R± (for the various angular modes l,m). 
These variables satisfy the decoupled equations 



d 2 R% n 
dr 2 



d 2 R l l n 
dt 2 



(3.13) 



with the effective potential 



Vl(r) = F 



1(1 + 1) -3M ±C + 4q 4 /i 



(3.14) 



where C — [9M 2 +4q 2 (l — l)(Z + 2)] 1 ' 2 . In turn, the gravitational and electromagnetic perturbations may be recovered 
from the fields R± (by certain linear combinations of the latter fields and their derivatives) [2D] . 

The potentials V± (r) are again cx F and therefore vanish at the two horizons, hence near the CH each of the fields 



R± l is decomposed into free ingoing and outgoing components, i.e. P u (u) + P v (v), as in Eq. fl3.I0|. Furthermore, the 
leading-order behavior of these two components in the early-CH domain u, v M was found [22 to be of the same 
form as in the scalar case, i.e. Eq. (3.11). As before, our objective is the proper-time variation of the perturbation 
fields, as recorded by a late-infall observer. The analysis of the previous subsections apply here with no modifications, 
implying that the variables R± undergo finite variations within proper-time intervals oc e~ K,Vch . The gravitational 
and electromagnetic perturbations (constructed from the variables i?^ n and their derivatives) are likely to yield a 
similar structure of an effective shock wave. 



3. Linear perturbations in Kerr spacetime 



In this subsection we very briefly address the case of a linearly-perturbed Kerr BH. The latter's internal structure 
is known to be similar in many respects to that of a RN BH. In particular there are two horizons, an event horizon 
located at r = r + and an inner horizon at r = r_ , where r± = M ± \J M 2 — a 2 . Not surprisingly, we find that the 
effective shock-wave phenomenon takes place in the Kerr case as well. 

We focus here on the gravitational perturbation, which is apparently the perturbation field of greatest physical 
relevance here. (Note that the Kerr background — unlike the electrovac RN background — admits pure vacuum gravi- 
tational perturbations.) 

The behavior of late-time gravitational perturbations inside a Kerr BH has been analyzed using two different 
formulations: (i) Analysis of metric perturbations (MP) h a p = g a/ 3 — g^g [9], (ii) Analysis pHHH] of the evolution 
of the Tcukolsky variables |23j tpo and "04. Both analyses examined the late-time gravitational perturbations, employing 
the so-called late-time expansion [T0J [TSl [24] . 10 They both focussed on the near-CH behavior, and led to similar (and 
mutually-consistent) results. 

For the sake of the present analysis the key result may be summarized as follows: Near the CH (that is, v — u 3> M), 
the linear MP decouple into a superposition of outgoing and ingoing components, namely 11 

h a p^Kp(u) + hlp(v). (3.15) 

Furthermore, h^a decays at v ^> M as an inverse power of v, hence for a late-infall observer h a p = h^Au) near the 
CH. 



Ref. [9] also considered nonlinear metric perturbations, that is, higher-order terms in the nonlinear perturbation expansion (which 
turned out, however, to be negligible compared to the linear metric perturbations). Ref. I10| also considered linear electromagnetic 
perturbations. In this section, however, we only consider linear gravitational perturbations. 

1 In the Kerr case we still define v = t + r» and u = t— r*, with r* now defined through dr/dr, = (r 2 — 2Mr + a 2 )/(r 2 + a 2 ) (where t 
and r are the Boyer-Lindquist time and radial coordinates). 
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In the Kerr background (unlike the RN case), infalling timelike geodesies may intersect the inner horizon r = r_ 
either at its ingoing or outgoing section. Here we shall focus on those orbits intersecting the outgoing section (these 
include, for example, all infalling geodesies with positive E and aL < 0). The effective shock- wave phenomenon will 
only occur for this class of orbits. 

Consider now an observer which falls into the Kerr BH (and heads towards the outgoing section of the inner 
horizon), at the late-infall limit. As before, we focus on the observer's history while moving between u — u\ and 
u = u 2 . Just like in the RN case, both proper-times T\ 2 scale at the late-infall limit as oc e~ KV ' h — and so does their 
difference At\2- Here k = n(M,a) is a certain positive constant, the inner-horizon's surface gravity. Again, we find a 
finite jump A/i Q ^ = h™p(v,2) — h^p(u\) in the metric, within an effectively-vanishing proper time oc e~ KVeh — namely, 
an effective gravitational shock wave. 

On the other hand, geodesies with aL > 2£Mr_ will intersect the ingoing section of r = r_. In the late-infall limit, 
these observers will hit the CH at its past boundary. 



C. Interpretation in terms of late-time Eddington frames 

We shall provide here a simple interpretation of the effective shock-wave phenomenon derived above. For con- 
creteness and simplicitly we present the explicit argument only for the RN case, but it applies to the Kerr case as 
well. 



The line element (2.6) preserves its form under a coordinate transformation of the form 

u — > ii = u — 5 , v — > v = v — <5 , (3.16) 

where S is any constant (this invariance reflects the time-translation symmetry of RN). We shall refer to different sets 
of Eddington coordinates — corresponding to different choices of 6 — as different Eddington frames (this terminology 
is borrowed from the analogous notion of Lorentz frames in Minkowski spacetime). Note that all tensors constructed 
from the metric are unaffected by this transformation. In addition, r* preserves its functional form, = (v — u)/2. 

Consider two infalling observers, which move along two identical worldlines related to each other by a time trans- 
lation. These observers cross the event horizon (EH) at Eddington times v\ h and v^ h , respectively (with > v^). 
Owing to this difference in v e h, the two observers will not share the same function v(t) (or u(t)). To brid ge thi s 



difference, we equip the second observer with its own Eddington frame (u,v), setting 5 — rffl — x& m Eq. (3.161 



"eh u eh 

.-,(2) - .,(1) 



Since now v^J = v^J, it is not difficult to show that v(t) of the second observer will be the same function as v(t) of 
the first observer — and the same relation will apply between u(r) and u(t). 



Consider now some linear perturbation field ty(u,v) on the background spacetime (2.6). Like all other tensorial 
quantities, \& is invariant under shifts in the Eddington frame. We assume that near the CH \t decouples to ingoing 
and outgoing components (like all linear fields considered above); and we shall be concerned here with the field's 
outgoing component, which we denote ^ u (u). 

Pick two u values 1*1,2, with the only requirement that A^f u = ^ 2 — ^1 7^ 0, where ^12 = ^ u (u — ^1,2)- We shall 
now examine how the second observer will experience this variation in from \&i to f 2, as a function of its own 
proper time. To this end we re-formulate the problem in terms of u rather than u. The change from ^1 to ^2 thus 
occurs while the second observer moves from u\ to U2, where 

- x r 1 (2) 

Let us now fix v^, and yet consider the late-infall limit for the second observer: — > 00. Evidently, in this 
limit both U\ and U2 are pushed toward 00. The corresponding (second-observer) proper times t(ui^) will thus be 
pushed to the same (finite) limiting value t(u —> 00), that is, the moment of inner-horizon crossing. In particular, the 
proper-time difference t(u2) — t(ui) will vanish in this late-infall limit. 

We conclude that at the late-infall limit, the finite variation in the perturbation field "J (which takes place 

between a certain pair of u values 1(1,2) occurs within a vanishing proper-time interval — at a moment which (at the 
limit) coincides with that of inner-horizon crossing. Thus we recover the effective Shockwave phenomenon for late-infall 
observers. 



IV. SPHERICALLY SYMMETRIC NON-LINEAR PERTURBATIONS 



In this section we shift our focus from linear perturbations on a RN (or Kerr) background, to nonlinearly-perturbed 
BHs. The main new ingredient is that now the infall orbit is disturbed by the MP, which in turn may influence the 
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observer's experience of the perturbation. For simplicity, we shall restrict attention here to a spherically-symmetric 
model of a nonlinearly-perturbed charged BH. We shall first present the model and describe the perturbed BH 
geometry, and then analyze the experience of late-infall observers in such a spacetime — demonstrating that the 
effective shock-wave phenomenon occurs in nonlinearly-perturbed BHs as well. 



A. Self-gravitating scalar field perturbations of a charged BH 



Let us consider a spherical charged BH perturbed by a spherically-symmetric self-gravitating scalar field. This 
model was investigated by several authors, primarily numerically [6HS] but also analytically [25J [26] (assisted by 
insights gained from earlier analytical investigations of the mass-inflation model [3]). The model consists of a 
masslcss, minimally-coupled, scalar field <j>, satisfying the covariant wave equation g a ^<j);ai3 — on the (self consistently- 
perturbed) metric of a spherically-symmetric charged BH. The scalar-field energy-momentum tensor 

Ta/3 = ~ ($;a<l>;P - ~0a/90;/t0 ! '*^ (4.1) 

acts as a source term in the Einstein equations (in addition to the electromagnetic contribution to energy-momentum) , 



yielding a system of nonlinear field equations for the metric functions [e.g. r(u,v) and f(u,v) in Eq. (4.2) below]. 
As for initial conditions, we consider here initial configurations wherein <f> is initially compactly supported outside 
the BH, as in 8 . (Alternatively, one may prescribe the initial data for (f>, including its presumed inverse-power tails, 
directly on the EH, as done in [7].) Evolving the initial data one then finds — not surprisingly — that at late times 
perturbations die out, and the BH settles down asymptotically to a member of the RN family, with charge q and a 
certain final mass M. The scalar perturbations decay as inverse-power tails. In particular, along the EH, <j> cx v~ n 
(typically with n = 3) at late times [HI [HUE]. These radiation tails fall into the BH and perturb its internal geometry. 

The perturbed metric in the BH interior is conveniently expressed in double-null coordinates. In particular, in 
Eddington-like coordinates 12 we write the line element in the form 

ds 2 = -f{u,v)dudv + T 2 (u,v)dn 2 2 . (4.2) 

The infalling scalar-field tail triggers the formation of a curvature singularity at the CH. This is a direct consequence 
of the infinite blue-shift that takes place at the inner horizon |27j . which leads to (almost) exponential divergence of 
the gradient of <f> — and of curvature. It turns out, however, that this is actually a weak [281 129] curvature singularity, 
located exactly at the CH (v — > oo). The metric tensor (in appropriate coordinates) approaches a finite, non-singular, 
limit as v — > oo. For example, we may use the Kruskal-likc coordinates 



with the line clement 



U=-e KU , V = -e~ KV , (4.3) 



ds 2 = H(U, V)dUdV + r 2 (U, V)dil 2 2 (4.4) 



(with H = KT 2 fe v ~ u \ note that U and V are both future-directed, and correspondingly H < 0). The CH is located 
at V — 0. Both r and H have finite, nonvanishing values at the CH. Yet, dr/dV diverges at V — 0, implying the 
presence of a null curvature singularity there (though a weak one). The scalar field (f> behaves in a manner similar to 
r: It is finite at V — > 0, yet d(j>/dV diverges at that limit. 

Perturbation theory predicts [25] (and numerical simulations [8j confirm) that in the early portion (i.e. u ^> M) of 
the CH, the metric functions deviate only slightly from their respective values in the unperturbed RN solution. The 
domain u,v 3> M is amendable to perturbative treatment. Correspondingly we express r and H as 

r(u, v) — r^jv(u, v) + Sr(u, v) , H(u,v) — H RN (u,v) + SH(u,v) , (4.5) 



The Eddington-like coordinate v may naturally be defined in the perturbed spacetime by using characteristic initial-value formulation, 
and setting r(v) and f(v) along the outgoing initial ray (beyond the end of the compact initial support of <f>) to be the same functions 
as in the unperturbed RN spacetime with the Eddington metric | |2.6| (setting F — > f therein). The key property of v is that it diverges 
at future null infinity, with / — > 1 on approaching the latter. An analogous procedure may in principle be applied to define u, yielding 
an ingoing null coordinate which diverges at the EH. 
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where the suffix "RN" denotes the corresponding function in the unperturbed RN spacetime. The perturbations 
<5r, 5H vanish in the limit u, v oo. (This limit corresponds to the past boundary of the CH, but also to t — > oo 
along spatial lines of constant v — u.) In the domain «,»» M, the scalar field is dominated by its linear perturbation 
term, and the MP Sr,SH by the second-order perturbation, as described below (see Appendix). 

As was mentioned above, r (like H and (f>) is finite along the CH. It initially starts at u — > oo with r = trn = r_, 
but subsequently shrinks steadily with time (— u), due to the focussing induced by the outflux of scalar-field energy- 
momentum across the CH. At some stage r shrinks to zero — at which point the null weak CH singularity terminates, 
and connects to a strong, spacelike, r = singularity [SHE]- 



B. Late-infall orbits 



We turn now to investigate the experience of late-infall observers in this spacetime of a nonlincarly-perturbed 
spherical charged BH. We consider (equatorial) infalling geodesies which are not necessarily radial. The angular- 
momentum parameter L is conserved in these geodesies, though E is no longer conserved. 

For concreteness we shall focus here on the behavior of the metric function r, which the infalling observer probes 
as a function of his proper time r. Physically, a rapid change in r will imply (for a finite-size observer) a rapid 
deformation in the tangential directions 9, ip. Again, we choose two u — const hypersurfaces, denoted u — ui 2 (with 
i*i > u 2 ), requiring that both hypersurfaces intersect the CH {before r shrinks to zero). While the observer progresses 
from u — u\ to u — u%, r changes by the amount Ar = r 2 — T\. Hereafter a quantity with a sub-index "1" or "2" will 
denote the value of this quantity as the worldline crosses the hypersurface u\ or ita, respectively. Since along the CH 
r is steadily shrinking, one finds that Ar =/= 0. 

Let us now evaluate the proper-time interval Ar = r 2 — t\ , using 



dx<*dx?\ 1/2 
dr=[-g a0 — — dU = 



\H\ 



,dV 
dU 



dip 
dU 



1/2 



dU. 



Note that H (like r) is bounded throughout the domain u\ > u > 112, so we can easily bound Ar by 



Ar < (H max ) x / 2 / 

JUi 



dV 
dU 



1/2 



dU, 



where H max denotes the maximal value of \H\ in the worldline's section between u\ and u^- The integral on the 
RHS is nothing but the proper-time length of a timelike curve connecting the two points (C/i, Vi) and (£/2,V2), in a 
fiducial two-dimensional spacetime with the flat metric ds 2 = —dU dV . It is bounded above by the (timelike) geodesic 

1/2 

connecting these edge points, whose length is (AV AU) ' , where AU = Ui — Ux and AV = V2 — V\ = e~ KVl — e~ KV2 . 



Clearly, AV < e 



< e 



1 , therefore 



Ar < (H max AU) 



1/2 e -M> ch /2 



Consider now the late-infall limit, which is the limit of large u e /v In this limit V\ : 2 — > 0. Therefore H max approaches 
r * all which is the maximal value of 
late-infall observers) in its final form: 



Hmaxi which is the maximal value of \H \ along the section u\ > u > u 2 of the CH. We obtain our bound on Ar (for 



Ar < Be- KV ' h/2 , 



(4.6) 



where B = (H^^AU) 1 / 2 is a parameter which depends on u\ and m 2 but not on the orbit's infall time. 

Wc conclude that late-infall observers will measure a non- vanishing variation Ar in the metric function r, within 
a short proper-time difference Ar which shrinks exponentially in the infall time v e h — which is again an effective 
shock-wave phenomenon. 



C. Do the late-infall orbits cross the CH? 



Our analysis so far did not make use of the perturbative nature of the metric field (at the early portion of the CH) . 
We merely assumed that the CH singularity is null and weak — and more specifically, that H admits a finite limiting 
value along the CH. Correspondingly, there was no need to restrict Ui and m 2 to the perturbative domain [u 3> M): 
We only required that the surfaces u = u\ and u = u 2 intersect the CH, rather than the spacelike singularity. However, 
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there still was one hidden assumption: We implicitly assumed that the late-infall observers will make it all the way 
from u\ to U2 without intersecting the CH (that is, with finite v). Once an observer intersects the CH, we cannot 
make any concrete statement about his subsequent experience, because the CH is by definition the boundary of the 
domain of unique prediction (for e.g. the metric functions). 13 

We therefore still need to complete this missing piece of the analysis. We shall show that as long as 1*2 is located 
in the weakly-perturbed domain of the CH (that is, U2 is large compared to M), the late-infall orbits indeed arrive at 
u = u 2 with finite v. 

The control on the growth of v will be achieved by monitoring the evolution of the geodesic's "energy" parameter 

E = -(u u + u v ). (4.7) 

Note that E (unlike L) is no longer conserved, because the perturbations destroy the exact ^translation invariance 
of the RN background. Yet, following the evolution of E will enable us to control u v , and thereby the evolution of v 
along the orbit. 



Equation of motion for E 



The lower-index covariant geodesic equation, applied to the Eddington-like metric (4.2 1, reads 

u„ = {l/2)g a ^u a u p = -(1/2) f^u u u v + (L 2 /r 3 )r. M . 

Therefore 

E = -(u u + u v ) = (l/2)(/,„ + /,,)uV - {L 2 /r 3 )(r >u + r, v ) . 

To get rid of the term u u u v in the RHS, we use the normalization condition g a pu a u^ = —1, which for the metric 
p~2| yields 

u u u" = (l + L 2 /r 2 )/f. (4.8) 

The quantity (f iU + f :V )u u u v then becomes 

(l + L 2 /r 2 ) [(In |/|), u + (In l/l),,]. 

Noting that the last term in squared brackets is equal to (ln|i/|) u + (In |-H]),u, we re-express E in terms of the 
Kruskal-like metric function H : 

E = - (l + ^) [(In |ff|),„ + (In \H\) tV ] ~(r „ + r, v ). (4.9) 

Note that in the background RN metric, both H and r are functions of v — u, hence E is conserved. Evolution of 
E will thus only result from the MP 5r and 8H . 

So far all equations were exact. To proceed beyond this point, we restrict attention to the early portion of the CH, 
where perturbations are presumably small, and carry out the analysis at the leading order in the MP. 14 Also, since 
we are dealing with the orbit's evolution very close to the CH, we may replace the background's functions and 
Hun by their corresponding inner-horizon values r_ and H_ = Hm\[(r = r_). We obtain 

E = A{8H :U + SH tV ) + B{5r, u + Sr tV ) = W(u, v), (4.10) 

where A and B are constants, 

1 ( L 2 \ L 2 



2H_ 

Recall that H_ (like r_) is a finite, non- vanishing constant. 



Furthermore, owing to the divergence of curvature at the CH singularity, it is unclear whether a classical extension beyond the CH will 
be physically meaningful. 

That is, we expand the various prescribed background functions (i.e. the functions of u, v obtained from the metric functions etc.) to 
first order in the MP. However, we do not linearize the worldline-related quantities, like E, E, u v , u u , etc. 
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2. Analyzing the evolution of E 



We re-writc Eq. (4.10) as 



dE/du = W{u,v)/u* 



(4.11) 



in which we view u w as a (yet unknown) function of the parameter u along the geodesic. One might choose to 
approximate this function u u (u) by the corresponding function for geodesies in the unperturbed RN geometry. We 
shall not proceed here in this way, because we do not want to assume a-priori that the accumulating perturbations in 
the orbital parameters must be small. Instead, we shall proceed by expressing u u in terms of E. To this end we use 
the contravariant version of Eq. (4.7), namely 



u u + u v = 2E/f. 



This, together with Eq. (4.8), constitutes a closed algebraic system for the two unknowns u u ,u v . One can of course 
write down the exact solution of this algebraic system. However, it will be more instructive to employ here the 
approximate solution, associated with the smallness of /: We are dealing here with the near-r_ region, where 15 
/ oc e K ^ a ~ v ' <C 1. The algebraic system thus yields the simple approximate solution (to leading order in the small 
parameter /) 16 



2E/f, 



(l+L 2 /r 2 )/(2E). 



Substituting this expression for u u in Eq. (4.11) we find 

d(E 2 )/du = Wf= [A(5H tU + 6H V ) + B(Sr, u + 5r. t 
Setting / = F = const x e K ( ll ~ t ') (where, recall, F = Jrn), we obtain 

d(E 2 ) 



du 



A{6H, U + 5H V ) + B(5r^ u + 5r, v ) 



k(u-v) 



(4.12) 



(4.13) 



(4.14) 



where we have absorbed the above const (in F) in A and B. Note that the RHS in this equation is a prescribed 
function of u and v (with no reference whatsoever to four- velocity) . 

At this point it will be useful to refer to the concrete form of the MP 5r and 5H . We focus here on the early portion 
of the CH (i.e. u 3> M), where the MP are small and decay as inverse powers of u and/or v. We denote by AE the 
modification in E acquired in the near-CH region, up to some it 2 3> M, due to the presence of MP. Based on the 
inverse-power form of the MP, in the Appendix we derive the bound 

\AE\ < C K/0 _(2 " +1) - C (v eh )-\ (4.15) 
where C is a certain parameter (independent of v e h). In particular, we find that for late-infall orbits |A2?| <C E. 



3. Analyzing the evolution of v 



We proceed now to analyze the evolution of v, showing that it remains finite throughout u > 112 (for any U2 3> M). 
From Eq. (4.12) we have 



dv _ u v 
du u v 



J_ 

AE 2 



(l + L 2 /r 2 ) 



(4.16) 



Now for late-infall orbits we already found that \AE\ <C E and hence we may regard E as constant (essentially the 
entrance value of E). Also we may set r = r_, and take the near-CH form of /, namely / = -e K ' u ~"' x const. We 
obtain 



dv 
du 



(4.17) 



15 Note that 6H <C -ffiyv (valid in the perturbative domain considered here) also implies that Sf (the perturbation in /) satisfies 5f <C fRN', 
that is, / S F oc e< u - v \ 

16 The algebraic system also admits a second solution , in w hich u u and u" are interchanged. However, since our late-infall observers enter 
the r ~ r_ region with large u u , it is the solution ( 4.12 I which actually takes place. 
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where C is some positive constant. Re-writing this as d(e KV ) / d{e KU ) = — C, we obtain 

e RV = -Ce KU + const . (4.18) 

[Notice that the last two equations are the same as those describing late-infall geodesies in the unperturbed RN 
geometry.] This expression is bounded above (by the const in the RHS). Thus, v remains finite throughout the range 
u 3> M — meaning that the orbit cannot cross the CH (located at v — > oo) in that domain. 



D. Concluding Remarks 



We found that the late-infall observers cannot cross the CH in the regime u ^> M . The analysis in section IV B 
then shows that the proper-time for these observers to move from u — u\ to u = 112, for any u\ > U2 S> M, decreases 
exponentially in the infall time v e h- Since r varies during this range by a finite amount Ar < 0, we inevitably face 
here the phenomenon of effective gravitational shock-wave: A discontinuity in the metric tensor, which propagates 
along a null hypcrsurface (the outgoing section of the inner horizon in our case). Physically, this means that any 
extended object will undergo a sudden deformation, by a certain amount, within an effectively-vanishing proper time. 

The amplitude of such a gravitational shock- wave may naturally be characterized by the (dimensionless) magnitude 
of the object's deformation. Specifically in the spherical model studied here, the shock's amplitude may be taken to 
be the dimensionless quantity |Ar|/r_. 

Two additional remarks are in order here: 



1. Since the discussion in section IV C assumed U2 S> M in order to treat perturbations perturbatively, one might 
mistakenly conclude that the amplitude of the gravitational shock wave can be weak (i.e. |Ar|/r_ <C 1). But 
this is not the case. It is clear from the above analysis that late-infall geodesies cannot fall across the CH until 



after the perturbations grow to be of order 1. We therefore see from section IV B that late freely- falling observers 



must face a strong gravitational shock wave whose amplitude is at least of order 1. 

2. In order to allow a simple discussion of the quantity E — ~{u u + u v ), which is conserved along geodesies in 
exact RN, we have so far assumed freely-falling worldlines. But let us now consider an accelerated late inf all- 
time observer. We choose the acceleration as a function of proper time to be bounded and such that, in the 
unperturbed RN geometry, the worldline would reach the outgoing inner horizon before crossing the ingoing 
inner horizon. Note that in the limit of large v e h all of these accelerations occur before reaching u = u\. As 
a result, in the unperturbed spacetime E becomes some constant Ef ina i for all u < u%. Corresponding late 
infall-time observers in the perturbed spacetime may thus be analyzed just as for the freely-falling observers 
discussed above but with the entrance value of E replaced by Efi na i. We conclude that any accelerated late 
infall-time observer who "would have reached the outgoing inner horizon in unperturbed RN" also experiences 
a Shockwave in the perturbed spacetime. 



V. EFFECTIVE SPACETIMES FOR LATE-INFALL OBSERVERS 



We have seen above that, in the limit of late infall times, observers who enter a perturbed Reissner-Nordstrom black 
hole experience an effectively unperturbed Reissner-Nordstrom geometry up to the point where they would expect to 
encounter an inner horizon at r = r_. At this point, observers who would have reached the outgoing inner horizon in 
exact RN then encounter a Shockwave across which the metric changes discontinuously. 

Describing the detailed nature of this discontinuity requires further investigation. While we will not attempt a 
precise treatment here, it is natural to expect that the above Shockwave in fact contains a curvature singularity, as we 
know that our observers will eventually reach such a singularity and we expect that, since they are already "nearly 
null" in the region described in section |PV] all proper times along their worldlines will be compressed to zero in the 
large v e h limit 17 . 

In addition, it is clear that some accelerated observers will reach the ingoing weak null singularity shown in figure 
[l] (left). For observers who enter the black hole very late, in the region where perturbations are very small, this may 
be accomplished using roughly the same set of accelerations (as defined in their own reference frame) as would be 



Note that this argument is equally valid if, for some kinds of matter, the spacetimes contain no spacelike singularity. In such cases the 
ingoing worldline will reach the weak null singularity, which is also a curvature singularity. 
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required to reach the ingoing part of the inner horizon in an unperturbed RN black hole. Furthermore, in the limit 
of late infall times, such observers arrive at the portion of the weak null singularity close to i + in figure [l] where the 
singularity is extremely weak, so that again such observers will measure no noticeable departure from unperturbed 
RN until r is exceedingly close to r_ . 

To compliment the precise calculations of sections III B and IV , we find it useful to give a brief heuristic argument 
that reinforces the above statements. To do so, note that taking our observer to enter the black hole at late times is 
equivalent to considering some fixed infall time and replacing the given spacetime with one in which the black hole 
formed (along with its full set of perturbations) at a much earlier time. As a particular example, we might consider a 
scenario which starts with empty Minkowski space in the interior and where both the matter that forms the black hole 
and our observer are dropped in from some large value of r at finite advanced times v ma tter and v b s respectively. For 
example, the observer might be ejected from a static space station. The matter could be dropped from a dense network 
of such space stations, or simply assumed to cross the r-coordiante associated with the observer's space station near 
the advanced time v ma tter- The observer is also given some particular instructions for executing accelerations as a 
function of proper time along his worldline and we take the matter to be released as some particular function of 
advanced time peaked around some v matter . Allowing v matter to vary leads to a one-parameter family of (isometric) 
spacetimes. We wish to study the limit v mat ter — oo holding v b s fixed. 

It is useful to assume that the perturbed spacetime preserves spherical symmetry. The advanced time v is thus a 
well-defined null coordinate everywhere in the globally hyperbolic part of the spacetime (shown at left in figure [I]) and 
the retarded time u is a similarly well-defined null coordinate in the region outside the event horizon. Furthermore, it 
is clear that two spacetimes differing by a shift of v ma tter can be described by essentially the same metric functions, 
albeit again with a corresponding shift of u, v. For example, introducing u = u — v ma tter and v = v — v ma tter, we may 
take the collapsing spacetimes to be described by the metric 



ds 2 = —f C oi{u, v) du dv + r 2 o[ (u, v)dCl 2 (5-1) 

for some fixed functions fcoi(u, v), r co i(u, v) which do not explicitly depend on v mat ter- Similarly, the center of spherical 
symmetry is described by some fixed curve u = c co i(v) and the event horizon is the curve u = +00. The fact that for 
any fixed v ma tter the spacetime approaches RN as u, v — > 00 implies that in the formal v ma tter ~°o limit with u. v 
fixed we have r col (u - v matte . r ,v - v ma tter) -> r(u, v) and f co i{u - v rnatter , v - v matt er) -> F(r(u, v)), where F{r) and 
r(u, v) are the functions defined in section [TT| for exact RN. 

We would like to find a similar construction for the region to the future of the event horizon. In that region the 
retarded time u cannot be defined simply by tracing outgoing radial null rays to / + . Rather than attempt to find a 
physically preferred definition of u, let us therefore consider any null coordinate u which is smooth inside the event 



horizon, has past-directed gradient V a it, and for which the metric again takes the form (5.1 1 on each member of our 
family of spacetimes where the functions f co i,r co i again approach F(r(u, v)) and r(u, v) as defined in section |ll| in the 
limit of large u,v with fixed v matt er- Note that inside the event horizon u will decrease to some w mQ tt er -independent 
minimum value u m i n at the point where the center of symmetry hits the singularity. We choose this u m i n to be finite. 
We expect that late-time observers experience an effective spacetime described by the formal v ma tter —00 limit 



(with u, v held fixed) of ( 5.1 1 for some definition of the internal u coordinate above 18 . Moreover, the essential features 
of this limit are independent of the particular definition of u. Indeed, by construction this limit is described by 
fcoi — > F(r(u,v)) and r co i — > r(u,v) on the domain u > u m i n + v maU er — > —00. All of the additional structure has 
been sent off to u — —00, which of course nevertheless lies at finite affine parameter along future-directed geodesies. 
Thus, as shown in figure [4j this limiting spacetime is just the region of unperturbed RN to the future of its past 
horizon and to the past of its outgoing inner horizon. However, rather than being smooth, the outgoing null surface 
u = —00 is a Shockwave into which all of the time-dependent structure of the spacetime has been compressed. We 
should also recall that, although it is no longer apparent from the limiting form of the metric, the entire ingoing 
surface v = +00 to the future of i + is a curvature singularity (i.e. the mass-inflation singularity) for any finite value 
of 



For example, if the metric coefficients of such a collapsing spacetime can be analytic in some null coordinates, then the action on the 
interior metric due to a shift of Vmatter will be determined by analytic extension of the action on the exterior metric, and this interior 
action can be used to define an interior null coordinate u that simply shifts with v mat ter- It is clear that the resulting u will have the 
properties required above near i + . 
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FIG. 4: Left: Some lines of constant u,v drawn on the spherical spacetime describing a collapsing shell perturbed by a 
massless scalar field. The long-dashed line (top right boundary) is the mass-inflation singularity. In contrast, the short dashes 
mark coordinate lines lying outside the physical spacetime (since r < 0). The vector field d u + d v is tangent to any surface 
u = ±oo or v = ±oo. The arrows indicate the direction of this vector field along such surfaces. The point marked i + is fixed 
under the action of this vector field. Right: The effective geometry for late-time observers obtained by flowing the interior 
of our spacetime backward along the vector field d u + d v , while similarly deforming the boundaries. As for figure 1 (right), it 
consists of the region of the unperturbed eternal Reissner-Nordstrom black hole with r > r_ , together with certain Shockwave 
singularities. The shock along the lower left boundary descirbes a jump from r = r+ to r = (corresponding to the regular 
center of spherical symmetry shown in the left diagram). Though it arises only in the formal late-infall limit, this shock is not 
actually accessible to any late-infall observer. In this sense the shock is fictitious. In contrast, the shock along the upper left 
boundary is accessible to late-time observers. It contains the final piece of the shell worldline and describes a jump from r = r_ 
to r — (corresponding to the spacelike singularity shown in the left diagram). 



VI. DISCUSSION 



Our work above examined the experiences of observers who enter spherically-symmetric charged asymptotically 
fiat black holes subject to both linear and non-linear perturbations. We also considered linear perturbations of Kerr 
black holes. Our emphasis was on the limit of late infall times, v e h — > oo. In this limit, our observers' observations 
agree precisely with those of similar observers in the unperturbed (stationary) black hole spacetime up to the point 
where the latter observers would reach the inner horizon. At that point, however, those late infall-time observers 
who would have reached the (smooth) outgoing inner horizon in unperturbed RN instead encounter an effective 
gravitational shock-wave in the perturbed spacetime. The shock's width (expressed in terms of the proper time of 
infalling observers) decreases with infall time as e~ KV " h , hence it effectively vanishes for late-infall observers. 

In such a gravitational shock-wave, the metric undergoes a discontinuity and changes by a finite amount 19 in 
vanishing time. An extended physical object hitting this Shockwave will thus undergo a sudden deformation (i.e. 
shear and/or contraction), by a finite amount, within an effectively-vanishing proper time. Owing to their short 
time scale, the gravitational tidal forces entailed in the shock will presumably dominate over all internal interactions. 
Thus, each single nucleon (say) will be deformed by a certain amount (typically of order 1) while traversing such a 
shock-wave. 

We also gave a heuristic argument that the experiences of late-time observers are described by the simple effective 
spacetime shown at right in figure |4j which in particular reproduces the above exact results. On the basis of this 
argument one expects the outgoing Shockwave to contain a strong curvature singularity 20 , with the area-radius r 



The magnitude of the shock-wave may also be infinite in some cases. In the discussion here we shall assume it is finite. 

For the case where the perturbations are associated with a massless scalar field. For more general matter fields, the Shockwave may 

contain only a weak curvature singularity. See footnote 1171 
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shrinking to zero across the Shockwave for the case described in figure [4| However, these final details remain to be 
verified by other techniques. 

While the explicit non-linear analysis in section |PV| was restricted to the spherically-symmetric case, we strongly 
expect that the same phenomenon of Shockwave formation will take place in nonsphcrical black holes as well, and 
particularly in perturbed spinning black holes. This expectation is based on the combination of several pieces of 
evidence: the experience of late-infall observers inside a linearly-perturbed Kerr BH (discussed in section IIIB), the 
sub-dominance of nonlinear perturbations in generically-perturbed spinning black holes (particularly near the early 
portion of the CH), and — more generally — on the profound similarity between the inner structures of charged and 
spinning black holes. 

It may be interesting to compare the properties of the two different types of null singularity that develop at the 
inner horizon: the shock-wave singularity at the latter's outgoing section, and the weak curvature singularity at the 
CH. We argue that from the physical point of view the shock-wave singularity is the more violent one — particularly 
for very late-infalling observers. To see this, consider two representative orbits in the RN geometry: Orbit (a) hits the 
ingoing section of r = r_ , whereas orbit (b) hits the latter's outgoing section. Let us now consider the experiences of 
these two observers on their approach to r = r_ , in case the RN spacetime is perturbed — focussing our attention on 
the late-infall orbits belonging to the two families (a,b). Observers of type (a) will hit a true curvature singularity at 
the CH. This singularity is weak, however, and any observer will only experience a finite integrated deformation up to 
the CH. Furthermore, this deformation will decrease with increasing v e h (as an inverse power thereof), and will vanish 
in the late-infall limit. On the other hand, observers of type (b) will experience, on approaching r = r_, a certain 
amount of deformation which does not decrease with increasing v e h- This is because the shock's amplitude approaches 
a certain limiting value at large v e h (it is only the shock's width which evolves at late time, it decreases exponentially 
in v e h). Thus, at least as far as the overall tidal deformation is considered, at the late-infall limit observers of type (a) 
will effectively feel no singularity at the CH, whereas late-infall observers of type (b) will feel a violent gravitational 
shock- wave singularity, whose magnitude is typically of order unity (at least). 

While our considerations were completely classical, our conclusions also suggest properties of quantum black hole 
microstates. In particular, the second law of thermodynamics implies that the late-time limit of any black hole 
spacetime should approximate the typical quantum state of the relevant ensemble. It is thus natural to conjecture 
that the typical state in an ensemble of charged or rotating black holes corresponds to the geometry displayed in 
figure [4] (right) together with quantum corrections, as opposed to say that shown in [4] (left). 

Some of our results carry over directly to generic (stable) black holes in any dimension. Indeed, in retrospect our 
linear analysis depended only on the general form of a (future-directed) time-translation dt between the inner and 
outer horizons of the unperturbed black hole, and in particular on the pattern of surfaces invariant under its action. 
The essential ingredients are just that: i) the (smooth) inner and outer horizons are invariant under <9 t , ii) the point 
^right m n g ure N (left) is an attractive fixed point of d t while i~l e t t is a repulsive fixed point, and iii) perturbations 
in the immediate vicinity of the outgoing inner-horizon decouple into independent u and v components 21 . Whenever 
these features arise, late infall-time observers who reach the outgoing inner horizon will experience any given linear 
perturbation as a Shockwave. 

At the non-linear level, the assumption that the perturbed spacetime approaches a stationary geometry at i^ ight 
which is an attractive fixed point of d t again implies that late infall-time observers experience an essentially un- 
perturbed solution up to the point where they would expect to reach an inner horizon. One is tempted to draw 
a conformal diagram similar to that in figure [4] (right) and again conjecture that infalling observers experience an 
effective outgoing Shockwave. However, investigating this conjecture in detail will require more sophisticated methods 
or, perhaps, numerical simulations. 

Throughout this work we have considered what one may call 'test observers,' which experience perturbations of 
the spacetime but do not source further perturbations. In contrast, any physical observer who falls into the black 
hole may be expected to create additional perturbations which now typically fall into the black hole at an advanced 
time comparable to v e h (and thus which cannot generally be ignored simply by taking the limit of large u e fc)- R 
would be interesting to understand what this implies for the experiences of such physical observers. For example, 
suppose our observer accelerates so as to reach the would-be ingoing inner horizon (the CH). Due to the decay of 
other perturbations as reviewed in section |IV A| for sufficiently large v e h the perturbations at the relevant part of the 
CH will generically be dominated by those sourced by our observer himself. One would expect him to encounter a null 
weak singularity whose strength (typically tiny, of order the mass ratio between the BH and observer) is determined 
by perturbations of his own creation. In contrast, even if our observer instead reaches the would-be outgoing inner 



At least for scalar fields with the usual kinetic term, one may derive (iii) from (i) and (ii) by defining V = XV and taking the limit 
A — > 00 to zoom in on the region near V = 0. Corresponding results may also hold for higher spin fields. 
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horizon, we see no reason for perturbations sourced by the observer to destroy the outgoing shock wave set up by 
earlier perturbations or to significantly change the experience described above. 

In our work above, we assumed that the spacetime approaches some non-extreme stationary black hole geometry 
near i + . To simplify the analysis, we assumed "generic" parameter values for this limiting BH, and in particular that 
it is not too close to extremality. While extending these results to nearly-extreme black holes is not difficult (see in 
particular footnote [24] in appendix A), it remains to understand the precisely extreme case. Taking the formal extreme 
limit of our late-time results (i.e., first taking v e h — > oo with r + 7^ r_ and then taking r_ — > r + ) agrees with the 
picture suggested in [3] (and supported by [3] ) in which late-time observers experience a singularity which effectively 
resides at the final event horizon. This suggests that such a picture should hold at least for horizons that approach 
extremality sufficiently slowly relative to the production of perturbations (and relative to the time scale set by the 
surface gravity). This in particular should apply for large black holes when the approach to extremality is due to 
quantum processes in cases where the extreme BH is quantum mechanically stable (e.g., with enough supersymmetry) . 
However, a complete analysis is again left for future work. 
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Appendix A: Bound on AE in slightly-perturbed RN 



In this Appendix we derive a bound on AE, based on Eq. ( |4. 14 1 and on the inverse-power form of the MP (at the 
early portion of the CH) . 

For a generic situation of a perturbed charged (or spinning) BH, the MP are dominated by the linear term, which 
is usually a superposition of a term cx v~ n and a term cx vT n , typically with n = 21 + 3 or possibly n = 21 + 2, 
where I is the multipolar number of the MP mode under consideration. However in our case the situation is different 
because we are only considering here spherically- symmetric perturbations, and there are no matter-free spherical MP. 
Therefore it is only the scalar field which is excited at the linear level. The MP are in turn sourced by the scalar-field 
energy- momentum tensor (4.1 1, which is quadratic in (derivatives of) <fi; hence the MP will first appear as (and will 



be dominated by) second-order perturbations. 

Thus, whereas in our case <f> is dominated by terms cx v~ n and cx u~~ n (with n = 21 + 3 = 3), the MP will 
be dominated by inverse-power terms of overall power —2n rather than n. Restated in other words, the MP will be 
dominated by a superposition of terms of the form CjkV^^u^ k , with non-negative integers j, k satisfying j + k = 2n = 6, 
where Cjk are certain constants. 22 



Correspondingly, the term in squared brackets in the RHS of Eq. (4.141 will be dominated by superposition of 
with non-negative integers f , k' satisfying j' 



terms c[-, k ,v 



k' 



2n + 1 (with some new constants c'-, kl ). 



Let us now restrict our attention to the contribution to d(E )/du coming from a single such term c'y k ,v 



-j',,-k' . 



d{E 2 ) [fkl] /du = dy k ,v" 



' U - k 'e K(v 



We can bound this contribution by 



-v) 



(Al) 



\d{E%, kl] /du\ < 



\Cj'k>\\ V n 



u~ k ' e KU , 



(A2) 



where u m m is the minimal value of v in the orbit's section under consideration. (u mi „ is always > v e h] Later we shall 
take v m i n to be the value of v when the orbit approaches the neighborhood of the CH.) 



To be more specific, assuming that <j> decays as v~ n along the EH, then perturbation analysis reveals [25] that in the neighborhood of 
the early portion of the CH Sr is dominated by a term oc v~ 2n ~ 1 plus another term oc u~ 2n_1 , and SH by a term oc v~ n u~ n (plus, 
possibly, additional terms oc v~ 2n and oc u~ 2n , though these two would be more sensitive to the choice of gauge). Note that as a 
consequence many of the coefficients </. /fc , (to be defined shortly) actually vanish, though some of them (for example c' n n+1 ) do not. 
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The term in the RHS of Eq. (A2) is an explicit function of u. Consider now its integral between a certain pair of u 
values Ui,uj satisfying Uj > Uf 3> M, We can now set v m i n to be Vi, namely, the value of v when the orbit approaches 
u = Uj. Also, since u^> M throughout the relevant domain, over the basic exponential time scale Su ~ 1/k the factor 
u~ k does not vary in an appreciable manner, therefore it may be pulled out of the integral. The latter is thus well 



approximated by just the integrand — namely the RHS of Eq. (A2) — divided by k (the error being smaller typically 
by a factor cx M/u). Thus, denoting the (j', k') contribution to the variation in E 2 (accumulated between Ui and Uf) 
by A(E 2 ) [j , k , ] , we find 



\A(E\, 



< 



IK)" 



u' k 'e KU 



Clearly this may be further bounded by the contribution at the upper limit (u = Ui), yielding the bound 



\A(E%, kl] \ < K - 1 |4 fc ,|(« 4 )^"(^) 



-k' ^{ui 



-Vi) 



(A3) 



The situation which concerns us is the motion of late-infall observers between two u = const hypersurfaces, u\ to 
it 2 (as described in Sections III|IV |. Our goal here is to show that the late-infall orbits will not cross the CH before 
approaching u — u 2 (for any u 2 ^> M). We therefore need to bound the variation in E 2 from the stage where the 
orbit approaches the neighborhood of the CH, up to u = u 2 . Correspondingly, in the above bounds on A(E 2 )j> k ' we 
should in principle set Uf — u 2 (though Uf actually drops out from our final bound), and choose Ui to be the value of 
u at which the orbit approaches the CH neighborhood. This notion of "approaching the CH neighborhood" is most 
naturally formulated by means of a certain value of r*, which we denote r*o- That is, a particular orbit is said to 
having arrived at the CH neighborhood at the point where = r*o, and we set our parameters u% and Vi to be the 
(it, v) values of that point (this in particular implies — Ui = 
value of at which e~ 2Kr * becomes <1. 23 

Substituting Ui — > Vi~ 2r*o in Eq. (A3), the exponent at the RHS reads e - 2fir *o. Also, since we are concerned here 
with the late-infall limit, in evaluating the factor (ui)~ k = (t'i — 2r*a)~ k we can safely ignore r*o compared to Vi 



2r* ). From Eq. (2.7), 7% is (roughly speaking) the 



(which is > u e /j), after which the RHS of Eq. (A3) becomes 



k V 



j'k' 



(v,)- { 3' +k,) e- 2Kr * 



(A4) 



Note, however, that the difference between Vi and v e h is, approximately, the quantity Av (defined earlier in Sec. 
Ill A), which is 0{M) (and independent of infall time), and hence is <S v e h in the late-infall limit. Therefore we can 
further replace (vi)^^ +k ' by (v e h)~ti +k K Also, since e _2Kr *° < 1 (it is in fact <C 1), we are allowed to entirely drop 
it from our bound. Recalling that j' + k' = In + 1, we obtain the bound on A(E 2s )yy k i^ in a much simpler form: 



\A(E%, 



< n- 1 \c' 



j'k' 



(Veh) 



-(2n+l) 



(A5) 



It remains to sum over the relevant pairs of integers j'k' (recall, these are certain non-negative integers satisfying 
j' + k! = 2n + 1). This summation is trivial and we get the bound on A(E 2 ) in its final form: 



|A(£ 2 )| <C'{v eh r (2n+1) ^C'{v eh )- 7 , 



(A6) 



where C = k 1 Y.j'\c'j, k , \ (setting k' = 2n + 1 — j'). Obviously the constant C is independent of v e h- 

We have thus established that at the late-infall limit, |A(i5 2 )| inevitably becomes <C E 2 . Equation (A6) thus also 
implies that AE (namely the modification in E accumulated in the near-CH region up to some ui 3> M) is bounded 
for late-infall orbits by 24 



\AE\ < C (v eh )-( 2n+ ^ =C(v eh )-\ 



(A7) 



Certainly there is an arbitrariness in the choice of r„o, but it does not affect our final results [e.g. Eqs. | |A6| and jA7[ ) below], which 
are independent of r »o- Note also that in the analysis here we do not attempt to include modifications in E caused by MP outside the 
BH, or even inside the BH but not in the CH neighborhood. These "non-CH" contributions do not concern us here, they can be shown 
by simpler arguments to decrease as an inverse power of infall time v e h. 

For simplicity, the above derivation assumed typical values of the black hole parameters such that k ~ 1/M; i.e., the black hole is not 
too close to extremity. We note in passing that a similar argument yields a bound of the same form for any re > 0. 
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where C= C'/(2E). 
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